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Abstract We study the magnetization curve of the Heisenberg model on the quasi- 
one-dimensional kagome-strip lattice that shares the same lattice structure in the 
inner part with the two-dimensional kagome lattice. Our numerical calculations 
based on the density matrix renormalization group method reveal that the system 
shows several magnetization plateaus between zero magnetization and the satu- 
rated one; we find the presence of the magnetic plateaus with the n/1 height of the 
saturation for n =1,2,3,4,5 and 6 in the S =1/2 case, whereas we detect only the 
magnetic plateaus of n =1,3,5 and 6 in the S =1 case. In the cases of n =2,4 and 6 
for the 5 = 1/2 system, the Oshikawa- Yamanaka- Affleck condition suggests the oc- 
currence of the translational symmetry breaking (TSB). We numerically confirm 
this non-trivial TSB in our results of local magnetizations. We have also found 
that the macroscopic jump appears near the saturation field irrespective of the spin 
amplitude as well as the two-dimensional kagome model. 
PACS numbers:75.10.Jm, 75.30.Kz, 75.45.+j 
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The kagome-lattice antiferromagnet has been extensively studied from both of 
experimental and theoretical approaches in the high expectation that the effect 
of the large quantum fluctuation and the strong frustration produce the exotic 
states. 1 1 2 i 3 1 4 i 5 1 6 i 7 1 8 i 9 1 10 1 11 From the viewpoint of the numerical studies, however, it 
is well known that the applications of the existing methods, the density matrix 
renormalization group (DMRG) 12 i 13 , quantum monte carlo (QMC) and exact di- 
agonalization (ED) methods, to the two-dimensional frustrated system such as the 
kagome antiferromagnet are difficult. The implementation of the DMRG method 
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Fig. 1 (Color online) Structures of the lattices: the kagome lattice (a), the quasi-one-dimensional 
kagome strip lattice (b).Antiferromagnetic bonds J (bold straight line) and ferromagnetic bond 
7f (dotted line). Sublattices in a unit cell of lattice (b) are represented by A, A', B, B', C, C, and 
D. 

is difficult in more than two dimensions although this method is very powerful 
for one- or quasi-one-dimensional systems under the open boundary condition 
and the QMC method comes across the so-called negative-sign problem for frus- 
trated systems. The ED method has the strong limitation of available system sizes 
although this method does not suffer from the limitation of the dimensionality 
nor negative sign problem. Due to these difficulties, there still remain many unre- 
solved problems about the nontrivial natures of the kagome antiferromagnets. In 
particular, very recent ED studies on the 5=1/2 antiferromagnetic Heisenberg 
model on the kagome lattice depicted in Fig. Qla) clarified that there exists an 
anomalous behavior at the 1/3 height of the saturation in the magnetization pro- 
cess, the magnetization ramp 910 whose critical behavior is quite different from 
conventional magnetization plateau and magnetization cusp. Therefore, the mag- 
netization curve of the kagome antiferromagnet is attracting the most attention. 

In this study, we investigate the magnetization curve of the Heisenberg model 
on the quasi-one-dimensional (Q1D) kagome strip lattice depicted in Fig.QJb) in- 
stead of the 2D lattice depicted in Fig.QJa). Note that the inner parts of the lattices 
in Fig. |TJb) are common to a part of the 2D lattice in Fig. QIa)- This model was 
originally introduced in ref. 14 for studying on the occurrence of the non-Lieb- 
Mattis ferrimagnetis m 15 ! 16 ! 17 in the ground-state of the 5=1/2 antiferromag- 
netic Heisenberg model on the spatially anisotropic kagome lattice 18 . However, 
the magnetization curve of this Q1D system has not been investigated in detail. 
We will show our numerical results obtained by the DMRG calculations not only 
in the S = 1/2 case but also S = 1 case. 

2 Model Hamiltonian 

The Hamiltonian of the present model is given by 

Jri? = ^^[Si,B • S,.C + S (j c • S/.D + S^c • Sj+^A + S^c • Sj+^B 

i 

+ S^c' ' S,-,B' + Sj C ' • S,-,a' + S^c' • Sm-i,d + S^c ■ S,- + i jB ' 

+ /S/,a • S/3 + S/,b • S;x> + S;x> ' S, )B ' + Sj b' ' Sj,A'] 
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i 

- fc£[S? A + .V; v + 5? B + .V; ,, + ,V a . + Sy + S? D ], (1) 

i 

where S,- ,s is an S = 1/2 or 5 = 1 spin operator at ^ -sublattice site in i-th unit 
cell. The positions of seven sublattices are denoted by A, A', B, B', C, C, and D 
in Fig.QJb). Note that the last term of eq.Q]is Zeeman term. The number of spin 
sites is denoted by N. Therefore, saturation magnetization value is M s = SN where 
M is equal to the z-component of the total spin S z tot . The number of unit cells is 
N/1; we consider N / 14 is an integer. Energies are measured in unit of J; we fixed 
/ = 1 hereafter. In what follows, we examine the magnetization curve in the case 
of 7f = — 1 by means of the DMRG method. 



3 Results 

We show our calculation results of the magnetization curves in Fig. |2ja) for S = 
1 /2 and in Fig. [2fb) for S = 1 . In the magnetization curve of S — 1 /2, we find 
6 plateaus with magnetization M/M s = 111, 211, 3/1, All, 5/7 and 6/7 where M s 
means saturated magnetization although we should pay careful attention to the 
finite-size effects especially under the open-boundary condition: in finite system, 
we regard the region of (= — ^) < M/M s < I as the 1/7 plateau, the region of 
§ < M/M s < + j|) as the 3/7 plateau, the region of f < M/M s < (f + f ) 

as the 5/7 plateau and the region of f < M/M s < (f + |) as the 6/7 plateau. 
In the magnetization curve of S = 1 , on the other hand, we find 4 plateaus with 
magnetization M/M s = 1/7, 3/7, 5/7 and 6/7. Note here that we regard the region 
of f < M/M s < (| + jj) as the 3/7 plateau, the region ofM/M s = (f + ^) as the 
5/7 plateau and the region of M/M s = + ^) as the 6/7 plateau in finite system. 
It is difficult to judge by using our present results whether or not the plateaus of 
M /M s = 2/1 and 4/7 exists in the thermodynamic limit. The issue of establishing 
the presences or absences of the 2/7 and 4/7 plateaus should be clarified in future 
studies. 

Oshikawa-Yamanaka- Affleck theorem 19 provides us the necessary condition 
for the magnetization plateaus as 

Q(S - m) = int. (2) 

where Q is the spatial periodicity of the wave function and m is the magnetiza- 
tion per site (m = M/N = MS/M S ). This condition tell us that the translational 
symmetry breaking should occur spontaneously in the magnetization plateaus of 
M/M s = n/1 where n is even numbe r 2Q i 21 . In order to confirm the symmetry 
breaking, we calculate the local magnetization (S^), where (A) denotes the ex- 
pectation value of the physical quantity A and c is the z-component of S, ^ . 

In Fig. O a ); we first present the correspondence relationships between each 
colored symbol and each sublattice ^ used in Figs. |3jb), [3jc) and [3jd). For ex- 
ample of the local magnetization on the plateaus without symmetry breaking, 
we show the calculation result of the local magnetization on the 1/7 plateau for 
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Fig. 2 (Color online) Magnetization curves of the kagome-strip lattice depicted in Fig. [2b). 
Panels (a) and (b) are the results of 5 = 1 /2 and 5=1 respectively. Double-headed arrows 
indicate the regions of each magnetization plateau in the N = 168 case for 5=1/2 and in the 
N = 84 case for 5 = 1 (see also the corresponding texts). 



the 5=1 case in Fig. ObO- One can immediately confirm the collinear spin- 
configuration reflecting the present situation that there are 7 sublattices in a unit 
cell of the lattice depicted in Fig. Ulb). In Fig. 3(c) and 3(d), on the other hand, 
we detect the oscillation in the local magnetizations within the same sublattices in 
the case of 2/7 and 4/7, respectively. The oscillation suggests that the spontaneous 
symmetry breakings occurs and that the wave function is degenerate. 

Finally, we discuss the relationships between our strip model and the original 
antiferromagnetic kagome model. We successfully find the macroscopic jumps 
near the saturation field in the magnetization curves of our strip models not only 
in the 5=1/2 but also in the 5 = 1 cases as well as in the case of the 5=1/2 
and 5=1 antiferromagnetic Heisenberg model on the original kagome lattice de- 
picted in Fig.QIa) 3 . This characteristic behavior was also observed in the case of 
some frustrated systems 22,23 and was proven in ref. 22. In our strip models, unfor- 
tunately, we does not observe clearly characteristic behavior of the magnetization 
ramp in the two-dimensional kagome-lattice antiferromagnet; Further examina- 
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Fig. 3 (Color online) (a)The correspondence relationships between each colored symbol and 
each sublattice | used in Figs. [3{b), [3jc) and[3hf). Panels (b), (c) and (d) are the calculation 
results of the local magnetization at each sublattice t, . 

tions are required to capture the behavior just outside the flat region with respect 
to magnetization. 

4 Conclusions 

We have studied the magnetization curves of S = 1/2 and S = 1 Heisenberg mod- 
els on the kagome strip lattice depicted in Fig.QJb) by the DMRG method. For S — 
1/2 case, we have confirmed 7 magnetic plateaus of M/M s =1/7,2/7,3/7,4/7,5/7 
and 6/7. In the case of M/M s =2/7, 4/7, 6/7, we confirm the occurrence of the 
translational symmetry breaking from our numerical results of the local magne- 
tizations. These symmetry breakings are suggested by the Oshikawa-Yamanaka- 
Affleck condition. For 5=1 case, on the other hand, we have confirmed 4 mag- 
netic plateaus of M/M s = 1 /7 ,3/7 ,5/7 and 6/7. We have also found the occur- 
rence of the macroscopic jump in these magnetization curves near the saturation 
field irrespective of the spin amplitude as well as the 5=1/2 and 5=1 anti- 
ferromagnetic Heisenberg model on the original kagome lattice depicted in Fig. 
Qla). 
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